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The values of renormalized Polyakov loops in the three lowest representations of 
5(7(3) were measured numerically on the lattice. We find that in magnitude, 
condensates respect the large-TV property of factorization. In several ways, the 
deconfining phase transition for N = 3 appears to be like that in the N = oo 
matrix model of Gross and Witten. Surprisingly, we find that the values of the 
renormalized triplet loop are described by an SU{3) matrix model, with an effective 
action dominated by the triplet loop. Future numerical simulations with a larger 
number of colors should be able to show whether or not the deconfining phase 
transition is close to the Gross- Witten point. 



1. Introduction 

't Hooft showed that the order parameter for deconfincmcnt in the SU{N) 
Yang-Mills theory is a global Z{N) spin given by a thermal Wilson line 
which wraps around the Euclidean time direction. The thermal Wilson 
loop in a representation TZ is given by 

\.n{x) = T exp (^g ' AI(x,t') dr'j . (1) 

Under a gauge transformation, 

LK(f) Vt^T^ix, 1/T) Lk(x) 0) . (2) 



1 
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There are important aperiodic gauge transformations: 

nN(^,l/T) = e''^ On(^,0) , 6'^-^ = !. (3) 

The Z{N) charge en (defined below) determines which loops can condense 
above the deconfinement temperature: 

(tr Lk) = , T<Td , enj^O (4) 

(tr Lk) ^ , T>Td , V . (5) 
The trace of the thermal Wilson line is the Polyakov loop ^, 

^n{x) = -^tvLnix) , (6) 
an 

and is gauge invariant. We divide by the dimension of the representation, 
d-jz, so that — > 1 as T ^ cx). 

The low-temperature, confined phase is Z{N) symmetric, so the expec- 
tation value of the fundamental loop vanishes below T^, while above T4 the 
Z{N) symmetry is broken spontaneously: 

(M = e'^ KMIt^O , e''^^ = l , T>Td. (7) 

The physics around Td is non-perturbative and it is necessary to employ 
lattice simulations. The expectation value of the Polyakov loop, however, is 
a bare quantity and contains ultraviolet divergences: the expectation value 
of the fundamental loop, and of loops in other representations, vanish in 
the continuum limit {Ni 00), at any fixed temperature T = l/(Nta). 
This is illustrated in Fig. 1. 

Renormalization of the expectation value of the bare Polyakov loop gives 
the exponential of a divergent mass times the length of the path 3,4,5,6,7,8. 

|(£N)|-exp(^-^j , (8) 

The rcnormalizcd loop, i-jz- is formed by dividing the bare loop by the 
appropriate renormalization constant, Z-ji- 

^n = ^iTi , ^7, = exp(-^) . (9) 

We have developed a method to extract the divergent masses non- 
perturbatively ^. The idea is to compute with a set of lattices, all at 
the same physical temperature, T, but with different values of the lattice 
spacing, a. The number of time steps, Nt = l/{aT), varies between these 
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Figure 1 . The bare fundamental loop as extracted from lattices with various Nt at three 
different temperatures. 



lattices and the divergent mass, am^"" , can be extracted by comparing the 
different values of the bare Polyakov loops. 

An alternate procedure was developed by Kaczmarek, Karsch, Pe- 
treczky, and Zantow ^ who obtain from the two point function of fun- 
damental loops at short distances. Our numerical values for the triplet 
Polyakov loop agree approximately with their values. So far, they have not 
considered loops in other representations. 

When the fundamental loop condenses it induces expectation values for 
all loops in higher representations, such as the octet and sextet loops. This 
can be precisely established in the limit of an infinite number of colors. 
Migdal and Makeenko observed that in SU (N) gauge theories, expectation 
values factorize at large N I0.ii.i2>i3,i4,i5,i6,i7^ infinite N, factorization 
fixes the expectation value of any Polyakov loop to be equal to powers of 
those for the fundamental and anti- fundamental, £^ = (^n)*, loop 



(10) 
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'N 



(11) 



This relation defines the Z[N) charge e-]^ of a loop in a given representa- 
tion. Thus, at infinite A'', any renormalized loop is an order parameter for 
deconfinement . 

To test these large N relationships numerically, for each loop we define 
the difference between the measured loop and its value in the large N limit. 
For three colors, the expectation value of the sextet difference loop is 

<54\ - (4) - (4)' , (12) 



(13) 



and that for the octet difference loop is 

2 



% =(4 - 



4 



Note that the difference loops vanish both at T — > and at T — > 00. If 
small, they indicate that factorization is approximately satisfied. 

2. Matrix Models 

Assuming that Wilson lines form the degrees of freedom in an effective 
theory, take as the partition function 

Z = j nidLN(i) exp (-<S(€k(*))) • (14) 

Here, i labels sites on a spatial lattice. The character expansion vastly 
reduces the possible couplings. Requiring the action to be Z{N) invariant 
allows a sum over Z{N) neutral loops: 

iV2 77^^7^(^) • (15) 

Add nearest neighbor interactions with couplings fin,'R.''- 

Sn = -— Si,^s;,^+f^'-° pn,n' Re (ndywii + n) • (16) 

Kogut, Snow, and Stone showed that for TV > 3 in mean field approxima- 
tion this model has a first order phase transition with a latent heat ~ 
and that 

KMI « ^ at T = T+. (17) 
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2.1. Three Color Example 

The simplest = 3 action includes just the triplet loop: 

S3 = -3/33 Re + ") • (18) 

Develop a mean field approximation by replacing all six nearest neighbors 
(in three space dimensions) by an average value, £0 = i^s)- 

Z = j dLg cxp(+18/334Re4) = exp(-9V) . (19) 
The mean field consistency condition is 

We have studied the condensates of the four lowest representations of 
SU{?>) as a function of the coupling (3^ ^, which are shown in Fig. 2. Our 




first observation is that the ordering (fundamental, octet, sextet) agrees 
with that obtained from the lattice. A fit of {(3){P3) obtained within the 
matrix model to {£3){T/Td) from the lattice gives 

= (0.46 ± 0.02) + (0.33 ± 0.02) ^ . (21) 
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Remarkably, the matrix model also predicts a sizable expectation value 
for the decuplet loop above T^; so far, this has not been studied on the 
lattice. Also note that expectation values of Z{3) neutral loops (octet and 
decuplet) are approximately zero in the confined phase. This is natural in 
a matrix model but does not follow automatically in an effective theory 
for a scalar £3 field 19^20,21,22,23 

While P3 is linear in the temperature, this doesn't seem to be true for 
all couplings. Fig. 3 illustrates this using the difference loops. Clearly, the 
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Figure 3. The difference loops computed from the matrix model and compared to the 
lattice data. 



"spikes" are much smaller and broader in the matrix model than in the lat- 
tice data. From this we conclude that for a more quantitative interpretation 
of the lattice results more terms are needed in the matrix model. 



3. Large— N Matrix Model and the Gross- Witten Point 

At infinite N, we do not have to consider loops in higher representations 
as independent degrees of freedom. In mean-field approximation, we have 
a single-site partition function: 

Z = / dLN exp (7V2(2/3£o) Re i^) = exp {-N^VawiPio)) , (22) 
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where (3 = (3n- The mean field condition is 

which amounts to minimizing the mean field potential 

Vmf{P,i)=Pi^ + VGwm . (24) 

This has been computed in the large TV limit by Gross and Witten 
The result is nonanalytic and is given by two different potentials: 

V-f=f}{l-f3)i' , (25) 
When P <1, the theory confines and £o = 0. For /? > 1 it deconfines and 



At /? = 1, the potential is completely flat, i.e. it vanishes identically for 
£e[0,i]. 

3.1. Large— N and Mass of Polyakov Loop 

The connected two point function of i is 

(£;,(x)£n(0)) - ^^^^^fe^ , IfHoo. (26) 

\x\ 

In the confined phase, m = cr/T, where a is the string tension. In the 
deconfined phase, one can define m oc niDebye- 
Computing vn? = d'^Vmf/d£^ gives 

m?. w 2(1 - /3) , /3 1" . (27) 

ml « V/9-1 , (3^1+ . (28) 
The string tension then vanishes at the transition as 

a(T) ~ (Td - T)i/2 , T^T^, (29) 
and the Debye mass, as 

mDebve{T)r.{T-Tdf/^ , T^T+. (30) 
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Hence, at the "Gross- Witten point" ^ there is a "critical" first order transi- 
tion whore the order parameter jumps at T^, yet both m_ and to+ vanish. 

Numerical results seem to indicate that the three-color Yang-Mills the- 
ory is close to the Gross- Witten point: i) the discontinuity of the renormal- 
ized fundamental loop is approximately 1/2; ii) the "spikes" of the difference 
loops are of order 1/A'' near and vanish at high T; iii) the string tension 
and the Debye mass drop sharply near the transition 

4. Conclusions 

It would be valuable to know from numerical simulations if the deconfining 

transition for more than three colors is close to the Gross- Witten point as 
well, or if that is unique to three colors. This is interesting, novel physics 
which can be obtained from lattice measurements at various N of the renor- 
malized Polyakov loops in hot SU{N) Yang-Mills theory. Further lattice 
simulations with higher accuracy will also constrain the couplings of the 
efi^ective matrix model description of the deconfining phase transition. 
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